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Abstract 

The conditions that allow us to consider the vacuum expectation value of the 
energy-momentum tensor as a statistical average, at some particular temperature, 
are given. When the mean value of created particles is stationary, a planckian distri- 
bution for the field modes is obtained. In the massless approximation, the temper- 
ature dependence is as that corresponding to a radiation dominated Friedmann-like 
model. 
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1. Introduction 

In the semiclassical approximation the vacuum expectation value (VEV) of the 
energy-momentum tensor (EMT) can be used as a source of the Einstein equations, 
in the form: 

Giiv = -8nG < T^ u >reg. (1) 

This approximation is useful as an asymptotic value of more complete theories 
where the gravitational field is also quantized, when the time is larger than the 
Planck time. The mean value of can be obtained using the expression: 

< out, 0\ f ^0, in > 
^ < out, 0\0, in > { ' 



which was introduced in ref.[l]. We can relate eq. (2) with the VEV with 
respect to some particular "in" or "out" vacua (see ref.[2]), i.e.: 

< >=< m,0|f^|0,m > +i ^ A *i T /^(0*m,i>0*m,j) 

=< out, 0|f M ,|0, out > +i Y^VijT^ ' out ^(j>* - outJ ) (3) 

where {<pi n } {J{4>*i n } and {4> ut} {J{4>* out} are the basis of solutions defined for 
the Cauchy surface labeled by "in" and "out" respectively. Ajj and Vij are func- 
tions of the Bogoliubov coefficients that determine the transformation between the 
two vacua. It is also useful to take into account that the expectation values can be 
expanded in terms of the vacuum polarization, which contains the local infinities, 
that are removed by the usual regularization methods (see ref.[3]) independently of 
the vacuum definition used , and a term which is function of the created particles, 
due to the interaction of the field with the curved geometry. This last term is of 
course null in flat space time. But in curved space it is in general infinite when it 
is calculated perturbatively. For instance in the particular case of fields minimally 
coupled to gravity, the infinity coming from the particle creation term, cannot be 
removed with the usual regularization methods for an arbitrary vacuum definition 
(see ref.[4]). However it is easy to see (ref. [5]) that when particle creation has a 
planckian distribution, which is lost when perturbative expansions are performed, 
the contribution to the EMT is finite and moreover is coincident with the standard 
cosmology which uses perfect fluid classical sources in the Einstein equations. A 
question that we can ask is; if it is reasonable to get a planckian distribution for the 
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mean value of created particles. The answer is afirmative because is obvious the fact 
that the actual background radiation of the Universe has a black body distribution 
with temperature approximately equal to 3°K (see ref. [6]). In another case, for 
example in the black holes studied by Hawking [7] , a planckian distribution is pre- 
dicted. Also the distribution seen by an accelerated observer (Rindler observer) is 
planckian like, therefore by the generalization of the equivalence principle we hope 
that for some Bogoliubov transformations we get a thermal spectrum. This result 
is shown in the present work for a scalar field minimally coupled to a Robertson- 
Walker metric. The expression obtained for the VEV of the EMT is compared with 
the one coming from the statistical average of T M „ at a given temperature, which is 
obtained using Thermo Field Dynamics (TFD) theory [9]. In a previous work [10] 
it was proved that a field in a curved geometry presents a form similar to that of 
the field in a thermal bath. In TFD the interaction with the bath is half-filled by 
the so called "tilde" modes, which is more natural that the interaction between the 
classical background and the quantum field. In the present work from, the compar- 
ison mentioned before, sufficient conditions on the Bogoliubov transformation are 
obtained, in order to get the same functional form for VEV of the EMT and the 
statistical average obtained from TFD. 



2. Time dependent Bogoliubov transformation 

In a pioneer work [11] Parker developed a consistent method to obtain the 

annhilation-creation operators at each time in curved space. The operators 
defined at different Cauchy surfaces are related by time dependent Bogoliubov 
transformations. Parker found in ref. [11] the differential equation satisfied by the 
Bogoliubov coefficients. This equation is other form to write the field equation. In 
order to resolve that equation it is necessary to give some initial conditions. The 
Bogoliubov coefficients give us the mean value of created particles between an initial 
time "to" an d a present time "£" . 

In order to make the paper more comprehensive, in this section we repeat in 
brief the formulation given by Parker in ref. [11]. 

The action for a minimally coupled scalar matter field is: 





and the metric is the spatially flat Robertson- Walker, given by 



ds 2 = dt 2 - a(t) 2 (dx 2 + dy 2 + dz 2 ) 



(5) 
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Therefore the variation of the action gives the field equation 

(V M ^ + m 2 )v? = (6) 

(with n = 0,1,2,3). 

The metric given by eq.(5) represents an unbounded space-time, therefore in a 
Fourier representation of ip integrals appear, in order to make the calculation easier, 
we can use a discretization in the same form as in ref.[ll], i.e. we can introduce 
the periodic boundary condition <p(x + nL,t) = <£>(x, t), where n is a vector with 
integer Cartesian components and L a longitude which goes to infinity at the end 
of the calculation (really in our case it is not necessary because the results used do 
not depend on L). Then we can introduce the set of functions {0k (x)} U {</>*k(;c)} 
defined by 

k (x) = — -— =expz(kx- / W{k, t')dt') (7) 
{La(t)y/ 2 V2W Jto 

with W an arbitrary real function of k =| k | and t. The field can be expanded 
in the form 

<p(x, t) = 5> k (*)0k(a;) + at k (^* k (a:)] (8) 

k 

where the operators a^(t) and a^(t) are time dependent and satisfy the com- 
mutation relations 

Mt),av(t)] =0, b\(t),a\,(t)]=0, [ak(t),a!v(t)]=6w (9) 
The operators act on the vacuum |0, t > in the form 

a k (£)|0,£>=0, a t k (t)|0,t>= \l k ,t > 
Moreover we define the operators 

^■k := a-k(t = t\) with t > t\ > to 

A k t :=a k t (t = ti) (10) 

Also we define the vacuum |0 >:= \0,t = t\ > on which the operators defined 
by eq.(10) act. These operators are related with the time dependent ones by the 
following Bogoliubov transformation 
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a k (£) = a(k,t)*A k + /3(k,t)A t _ k (11a) 

a_ k t(t) = a(k, £)^_ k + /3(k, t)*A k (116) 

clearly we use the conditions 

a(fe,*i) = l, /5(fc, ti) = (12) 

We can rewrite the field y?(x, t) in terms of the new operators in the form 

p(x, £) = ^[A k ^ k (x) + At k ** k (a;)] (13) 
k 

with \I/ k (x) and solutions of eq. (6). We can make the separation 

*k(aO = M k > t) exp zk.x (14) 

with 

, /, s 1 r „ n* -*f t Wdt' .* if t Wdt\ 

(La(t)) 3/2 (2W(k,t))V2 L 
Replacing eq. (14) in eq. (6) we have 

/ i f Wdt' . p -i f Wdt'x,^ W / . j T vKdt' . a -i f Wdt'\ 

ae J + pe J )M — — (ae J + pe J ) 

■ o-tja/ • i f Wdt' ib —i f Wdt'\ .-if Wdt' , n -i f Wdt' n /-, r\ 

+2iW (ae J — (3e J ) + ae J + pe J =0 (15) 

with 

M := + ^(W/W) 2 - \{a/af - \{a/a) + u 2 - W 2 (16) 

In the following we will call uj 2 = k 2 /a 2 + m 2 . 

From eqs (12) and (13), Parker [11] proves that the Bogoliubov coeficients 
satisfy the equation 

p=-aexp(2i[ W(k,t')dt') (17) 

J to 
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From the functional form (14), for which the invariance of the Klein-Gordon 
product holds, we obtain: 

\a\ 2 -\(3\ 2 = l (18) 

Using eq.(ll) we note that := |/3k| 2 =< 0\a\(t)ak(t)\0 > is the mean value 
of particles created between t\ and t. Eq. (18) allow us to write (as in ref.[ll]): 

a(k, t) = (M) cosh 0(k, t) (19a) 



P(k, t) = e*™^'*) sinh6>(k, t) (196) 

7a 5 7/3 and are functions of k and t that must satisfy the initial conditions 
given by eq.(12) and the field equation (15). As we can see from eq. (19 a) the 
number of particles is only a function of 9, i.e. ; = sinh 9. 

Replacing eqs (19) and (17) in eq. (15) and splitting in real and imaginary 
parts, the following system of equations is obtained: 

(1 + tanh#cosr)M + 2W7 Q = (20a) 



MsmT + 2W9 = (206) 
with r := 7a + 7/3 - 2 / t * Wdt' 

The solutions to eqs (20) give us the set of time dependent Bogoliubov trans- 
formations compatible with the Robertson-Walker metric. The stationary particle 
creation is a particular case with 9 = 0, therefore from eq. (20 b) we have: 

M = (21) 



3. VEV of EMT as a thermodynamical object 

In ref.[4] the analogy between the Bogoliubov transformation which relates 
states at different temperature with those at different Cauchy surfaces was shown. 
Here we will exploit this analogy in relation with the VEV of EMT. First of all we 
calculate the classical EMT in curved space-time in the usual form (see ref.[2]), i.e. 
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we obtain the expression 

T »v = \ <Pu} - ^9i*v{<P a , <Ar} + ^m 2 g flu {cp, if} (22) 

where { , } is the anticommutator and ip^ := d^tp. 

We can now calculate the EMT operator replacing eq.(13) in eq.(22) obtaining 



kk' 

+A\A v D liV [*\(x), Vv(x)} + A^A^,D^[y\(x),y\,(x)} 

+{k <-> k'}} (23) 

with {k <-> k'} indicating that the last term of eq.(23) is equal to the first 
changing k by k'. The differential operator is defined by 

D^[cp, if>] := -d^ipdvip - ^g lxl/ d a (pd a i() + -m 2 g IJiU ip'ip (24) 

We can now calculate the VEV of T^ v using the vacuum |0 >, and obtain: 

< 0|T^|0 >= 2^ k i?e{ZV[* k (x), V\(x)}} (25) 
Replacing now eqs (11) and (13) we also have 

^(x) = a*Mx)+P*<t>-* (26) 
Putting this last equation in eq. (25) yields 

0|f^|0 >= 2^i?e{(l + 2\(3\ 2 )D^[^ k *] + 2a*f3D„ v [^ 0_ k ] 



< 



k 

+ I(|/3| 2 ji^[«/> k , </> k *] + ^(a'flE^fa, 0_ k ] 
+\(S 0li S 0v - \g^)[(\a\ 2 + 0\ 2 )\^\ 2 + 2a*/?*0_ k k ]} (27) 

with 

if) (28) 
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For simplicity we can consider the case in which the particle creation is sta- 
tionary, then $ = a = 0, therefore we have 

< 0|f M ,|0 >= 2^e{(l + 2|/3| 2 )ZV[</> k , k *] + 2a*pD^[fa, 0_ k ]} (29) 

k 

As we will see the difference between the VEV given by eq. (29) and a statistical 
average is the "interference" term, for the modes k and — k, of eq. (29). 

A simple way to calculate the statistical mean value of some physical observable 
is by means of the thermo field dynamics (TFD) formulation (see ref. [12]). In this 
formulation the idea is to take into account the interaction between the system 
and the thermal bath by means of quantum fluctuation of the bath (see also ref. 
[9]), which are called "tilde" fields. These fields are auxilliary because they are 
not measurable. Those fields only act in the thermalization of the system modes. 
The tilde modes are associated to tilde creation annhilation operators which act on 
the vacuum of the thermal bath |0 >. The space of states, in this formulation is 
extended in order to include the tilde and the states of the system, i.e.: 

{\n, h >} = {\n >} <g> {\h >} 

We will call the total vacuum ||0 >:= |0 > <E>|0 >. Moreover we introduce a 
T parameter in order to identify the temperature of the state, i.e. |0,T > and a 
Bogoliubov transformation that relates the vacuum at zero temperature with the 
one at temperature T: 

a k (T) = A k cosh 9(k, T) — sinh #(k, T) (30) 

In order to obtain the zero temperature operators, #(k, T = 0) = is necessary. 
Then we have the creation annhilation bosonic fields A k and A^ k which operate 
in the form 

A k ||0 >= 0. 
^kHO >= ||l k >, etc. 
with the commutation relation 

[Ak, = 5 kjk / 

The tilde operators, represent the quantum effect of the reservoir and satisfy 

ik||0 >= 0. 
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AjjO >= ||l k >,etc. 
[A^Al] =5 k)k , 

The thermal operators a k (T), a^ k (T), a k (T), aj^T), which satisfy: 

a k (T)||0,T>=0 

a t k (T)||0,T>= ||l k ,T>,ete. 
[a k (T), a t k /(T)] = <5 kjk / 

and in analogous form the tilde operators. The tilde and no tilde operators are 
mutually commutative. 

With the new Fock space we can calculate the statistical mean value of any 
physical observable, for example if we have the physical operator A, the mean value 
can be obtained by means of: 

< A >=< 0,T||A||0,T> 

We will calculate now the statistical mean value of the energy momentum 
operator T^ v . In orther to do that we can represent the operator as a functional of 
the thermal modes k r an d 4>*\iT in a wa y analogous to that in curved space time: 

= J2i a *( T ) a v( T ) D ^[Mx),<i>VT(x)} + a ]i (T)a\ / (T)D llu [ ( p ]iT (x), ( p\, T (x)] 
kk' 

Wk(T)a k ,(T)^[0\ T (x),^ T (x)]+at k (T)at k ,(T)^[0\ T (x),0* k , T (x)] 

+{k^k'}} (31) 

It is easy to prove (see ref. [13]) using the transformation given by eq. (30) 
and the inverse that 

< 0\\a^(T)ak{T)\\0 >=< 0, T\\ A k U k ||0, T > 

and in analogous form with ak(T)ak(T) and ak(T)ayJ(T). Therefore we can do 

< T„ v >=< 0||T^[a k (T), a k f (T), a v (T), a^(T)] ||0 > 

then the following expression is obtained 

< 0|V|0 >= 2^Re{(l + 2n k )D lxu [ ( j) kT , kT *]} (32) 

k 
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where 

n k =<0||a k t(2> k (r)||0> 

As we can see from eq.(32) the 0k modes will be equivalent to the thermal ones if 
the interference term in eq. (29) is null. 

The difference with the "genuine thermality" has the same cause as that in 
the case of the Rindler vs Minkowski observers, analized in ref. [14]. It is related 
with the fact that < 0|ak(i)ak'(t)|0 >^ 0. In ref. [14] the thermality is restored 
assuming that the observation is restricted to a finite spacetime region, which is 
introduced mathematically by means of the wave packets (see also ref. [8]). In our 
case instead of using wave packets we will exploit the fact that the interference term 
is oscillant. First we calculate each component of that term: 

M2a*/3ZV</>k,^-k]} = ^cosh#sinh#|0| 2 {[^ + 3#) 2 

+u 2 - W 2 ]cosT- (^ + ?>H)Ws\nT} (33) 
i2e{2a*/?D O ;[0k,0-k]} = -Re^ M^kj} (34) 
i2e{2a*/3Djj[0k,0_ k ]} = -cosh0sinh0|0| 2 (2/c 2 ,,- - a 2 w 2 )cosT (35) 

There are not contributions coming from the term given by eq. (34), because when 
the summation X]k=-oo + °° * s performed, a self cancellation is produced. 

The terms given by eqs (33) and (35), when the limit to the continuum is taken 
(X)k — *■ (L/2ii) 3 J <i 3 k), are proportional to the integrals: 

7 o= I k 2 dk{[h^- + 3H) 2 + u 2 - W 2 ]cosr - (^ + 3#Wsinr}cosh#sinh# 
J AW W 

(36) 

poo 

Ijj = / k 2 dk{2k 2 - a 2 w 2 )cosrcosh#sinh# (37) 
Jo 

(where we used d s k = Ank 2 dk, and V as in eqs (20)). The particle creation is 
not modified when 7 changes, as we can see from eq. (19). Then we can propose a 
7 phase with the functional form 

7 = fj,k (38) 

If we take \x — > 00 and introduce an ultraviolet cutoff (which is less restritive than the 
wave packet and can be justified because k/a is the physical momentum associated 
to the k mode, as is shown in ref. [11]) in the integrals (36) and (37) we can then 
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apply the Riemann-Lebesgue theorem, which is as follows: "If / is a real function 
absolutly integrable in the interval [a, b], then \im. 1 ^ 00 j^f(t) sin(7t + 8)dt = with 
5 a real constant" . Clearly this is also true when we have "cos" instead of "sin" . 
By means of this theorem we can eliminate the interference term in eq. (29), and 
it results equivalent to the statistical mean value given by eq. (32). 

4. Thermal spectrum condition 

In order to obtain thermal distribution for the k modes, we need to apply 
an extremun condition to a thermodinamical potential. This potential will be a 
function of the energy and the entropy of the bosonic gas in curved space-time. 
The energy can be calculated by means of the metric hamiltonian H in the form 

E =< 0\H\0 > 

where 

H = J a 3 d 3 xT bo 

in the discretized formulation is J d 3 x = L 3 . For eq.(29) when the oscillatory term 
is dropped, we have: 

E = J2(\ + n ^ ( 39 ) 

k 

with 

fk =2l7 l I ( F + 3H)2 + l "' 2 + w2 ) (40) 
the energy by mode. As the energy for the system studied has the form of a set 
of decupled quantum oscilators, it can be considered as a bosonic gas in flat space 
time, therefore we can use the expression for the entropy used in Minkowski space. 
In order to do that we can write the operator (see ref.[12]): 

K = — ^{4(t)a fc (t) logsinh 2 6> - a k (t)al(t)cosh 2 9} (41) 
k 

The expectation value of the operator given by eq.(41) is the entropy of a bosonic 
gas [15]: 

K =< 0\K\0 >= - J2i n k l ogn k - (1 + n k ) log(l + n k )} (42) 

k 

In a way analogous to ref.[12] we will get an extremun, with respect to the Bogoli- 
ubov angle 6, of the thermodinamical potential 

A = —TK + E 
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where T is a Lagrangian multiplier that can be interpreted as the temperature. 
Then performing the variation 

5A/56 = 



the following eq. is obtained 



, sinh 2 9 

T log j- + e = 0. 

cosh 9 



therefore we get a planckian spectrum 

"k = -jj— (43) 
e t — l 

(whit units so that Ti = c = ks = 1) 

It is interesting to note that we can eliminate the interference term in the 
hamiltonian, without using the phase condition given by eq. (38), by means of the 
criterium known as hamiltonian diagonalization. As we can see from eq. (33) that 
condition is 

W 

w = ~ 3H (44a) 

W 2 = u 2 (446) 

This condition is, moreover, compatible with the dynamic equation for W, given by 
eq. (21). When at the observation time holds eq. (44), Doot'/'k, 0-k] = is assured. 
Of course the condition (44) has a meaning when the calculation is "in-out" , with 
the parameter "t" fixed at the observation time. Then we can think of condition 
(44) as Cauchy data on W and W. So when condition (44) is satisfied the energy 
by mode is (as we can see from eq. (40)). 

ek = u (45) 



In particular if we suppose that the rest energy can be dropped against the 
kinetic one, i. e. k 2 /a 2 » m 2 , we can do ek — k/a, and obtain 

^k = — !— - (46) 

e^T — 1 

From eq.(39) going to the continuum set of modes, the particle creation con- 
tribution to the energy density is given by 



2n 2 a 4 



DO 



1 < ,3 



k 6 n k dk (47) 
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Replacing eq.(46) in eq.(47) and performing the integral, we have 



P = 



2tt 2 



1 



rp4: 




dk' = —T 
30 



(48) 



which is the expression obtained from the standard phenomenological radiation 
dominated Friedmann cosmology [16]. 

5. Conclusions 

The VEV of the EMT is different to the statistical mean value, even for the 
stationary particle creation. In the last case the difference is an interference term 
between the modes k and — k, which can be eliminated by a condition on the phase 
of the Bogoliubov transformation. At the hamiltonian level the coincidence with 
the statistical value is satisfied by means of the criterion known as diagonalization 
of the hamiltonian condition (that for scalar field is identical to the condition of 
minimization of the energy ). By means of this condition also the mean value 
< T 00 >, defined by Utiyama-De Witt [1], and the VEV of the EMT are coincident 
(see eq.(3)). 

The coincidence between the VEV and the statistical mean value, is another 
view of the analogy between the Bogoliubov transformation that connects the dif- 
ferent vacua in curved space-time and the one that relates vacua to different tem- 
perature in TFD. 

From the extremal conditions on the thermodynamical potential, we obtained 
a thermal spectrum. The fact that particle creation turns to be thermal avoids the 
existence of infinities coming from this effect, in the source of the Einstein equations. 
This is true also when the hamiltonian diagonalization conditions are used. In 
general those conditions produce inconsistences with the adiabatic regularization 
[4], [17]. All of that suggested the utilization of the usual covariant regularization 
tecniques, in order to eliminate the infinities coming from the vacuum polarization 
term, together with some thermodynamic conditions that "regularized" the creation 
particle term. That approach allowed us to take into account the back-reaction effect 
in an easy way, without using numerical calculation (see ref.[5]). 
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